Casimir effect for curved geometries: PFA validity limits 
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We compute Casimir interaction energies for the sphere-plate and cylinder-plate configuration 
induced by scalar-field fluctuations with Dirichlet boundary conditions. Based on a high-precision 
calculation using worldline numerics, we quantitatively determine the validity bounds of the prox- 
imity force approximation (PFA) on which the comparison between all corresponding experiments 
and theory are based. We observe the quantitative failure of the PFA on the 1% level for a curvature 
parameter a/R > 0.00755. Even qualitatively, the PFA fails to predict reliably the correct sign of 
genuine Casimir curvature effects. We conclude that data analysis of future experiments aiming at 
a precision of 0.1% must no longer be based on the PFA. 
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Measurements of the Casimir force pj have reached a 
precision level of 1% USSIESSIi. Further improve- 
ments are currently aimed at with intense efforts, owing 
to the increasing relevance of these quantum forces for 
nano- and micro-scale mechanical systems; also, Casimir 
precision measurements play a major role in the search 
for new sub- millimeter forces, resul ting^ in important con- 
straints for new physics 0, El 0, ILlll3l| . 

On this level of precision, corrections owing to material 
properties, thermal fluctuations andgeometry dependen- 
cies have to be accounted for 0, 0, 0, . In order 
to reduce material corrections such as surface roughness 
and finite conductivity which are difficult to control with 
high precision, force measurements at larger surface sep- 
arations up to the micron range are intended. Though 
this implies stronger geometry dependence, this latter ef- 
fect is, in principle, under clean theoretical control, since 
it follows directly from quantum field theory [Tsl ]. 

Straightforward computations of geometry dependen- 
cies are conceptually complicated, since the relevant 
information is subtly encoded in the fluctuation spec- 
trum. Analytic solutions can usually be found only for 
highly symmetric geometries. This problem is particu- 
larly prominent, since current and future precision mea- 
surements predominantly rely on configurations involving 
curved surfaces, such as a sphere above a plate. As a gen- 
eral recipe, the proximity force approximation (PFA) |l9j 
has been the standard tool for estimating curvature ef- 
fects for non-planar geometries in all experiments so far. 
The fact that the PFA is uncontrolled with unknown va- 
lidity limits makes this approach highly problematic. 

Therefore, a technique is needed that facilitates 
Casimir computations from field-theoretic first princi- 
ples. For this purpose, worldline numerics has been de- 
veloped |20j . combining the string-inspired approach to 
quantum field theory [21j with Monte Carlo methods. As 
a main advantage, the worldline algorithm can be formu- 
lated for arbitrary geometries, resulting in a numerical es- 
timate of the exact answer [22| • For the sphere-plate and 
cylinder-plate configurations, also new analytic methods 
are currently developed and latest results including ex- 
act solutions are given in 0,13]. In either case, quanti- 
tatively accurate results for the experimentally relevant 



parameter ranges are missing so far. 

In this Letter, we use worldline numerics |2(],|22] to ex- 
amine the Casimir effect in a sphere-plate and cylinder- 
plate geometry for a fluctuating scalar field, obeying 
Dirichlet boundary conditions ("Dirichlet scalar"). Wc 
compute the Casimir interaction energies that give rise 
to forces between the rigid surfaces. Thereby, we quan- 
titatively determine validity bounds for the PFA. Apart 
from numerical discretization, for which a careful error 
management on the 0.1% level is performed, no quantum- 
field-theoretic approximation is needed. 

We emphasize that the Casimir energies for the Dirich- 
let scalar should not be taken as an estimate for those 
for the electromagnetic (EM) field, leaving especially the 
sphere-plate case as a pressing open problem. Never- 
theless, the validity constraints that we derive for the 
PFA hold independently of that, since the PFA approach 
makes no reference to the nature of the fluctuating field. 
If an experiment is performed outside the PFA validity 
ranges determined below, any comparison of the data 
with theory using the PFA has no firm basis. 
Casimir curvature effects. - An intriguing property of 
the Casimir effect has always been its geometry depen- 
dence. As long as the typical curvature radii Ri of the 
surfaces are large compared to the surface separation a, 
the PFA is assumed to provide for a good approximation. 
In this approach, the curved surfaces are viewed as a su- 
perposition of infinitesimal parallel plates fill ITgj . The 
Casimir interaction energy is obtained by an integration 
of the parallel-plate energy applied to the infinitesimal 
elements. Part of the curvature effect is introduced by 
the choice of a suitable integration measure which is gen- 
erally ambiguous, as discussed, e.g., in J2{|. For the case 
of a sphere with radius R at a (minimal) distance a from 
a plate, the PFA result at next-to-leading order reads 

E PFA (a, R) = E^ A (a, R) (l {j}-^ + 0((-|) 2 )) , (1) 



Ejft A (a,R) = -c PP 
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where the upper (lower) coefficient in braces holds for 
the so-called plate-based (sphere-based) PFA. They rep- 
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FIG. 1: Contour plot of the negative Casimir interaction en- 
ergy density for a sphere of radius R above an infinite plate; 
the sphere-plate separation a has been chosen as a — R here. 
The plot results from a pointwise evaluation of Eq. using 
worldlines with a common center of mass. 



resent two limiting cases of the PFA and have often been 
assumed to span the error bars for the true result. Fur- 
thermore, cpp = 2 for an EM field or a complex scalar, 
and cpp = 1 for real scalar field fluctuation. 

Heuristically, the PFA is in contradiction with Heiscn- 
berg's uncertainty principle, since the quantum fluctu- 
ations are assumed to probe the surfaces only locally 
at each infinitesimal clement. However, fluctuations are 
not localizablc, but at least probe the surface in a whole 
neighborhood. In this manner, the curvature information 
enters the fluctuation spectrum. This quantum mecha- 
nism is immediately visible in the worldline formulation 
of the Casimir problem. Therein, the sum over fluctu- 
ations is mapped onto a Feynman path integral. Each 
path (worldline) can be viewed ("IS i\ random spacetime 
trajectory of a quantum fluctuation. Owing to a generic 
spatial extent of the worldlines, the path integral directly 
samples the curvature properties of the surfaces ■ 

For the Dirichlet scalar, the worldline representation 
of the Casimir interaction energy boils down to (2^, H|| 



1 1 
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^e" m ° T (ex[x)) x . (3) 



The expectation value in (j^J) has to be taken with respect 
to an ensemble of closed worldlines, 



Vx 



. e 



(4) 



x(T)=x(0) 



with implicit normalization (l)^ = 1. In Eq. 6s[x] = 
1 if a worldline x intersects both surfaces £ = Si + £2, 
and 9s[s] = otherwise. A worldline with 9j[i] = 1 
represents a boundary-condition violating fluctuation. 
Its removal from the set of admissible fluctuations con- 
tributes to the negative Casimir interaction energy. 

We evaluate the worldline integral with Monte Carlo 
techniques, generating an ensemble of tt-l worldlines with 
the v loop algorithm [22| ■ Each worldline is characterized 
by N points after discretizing its propertime. In this 
work, we have used ensembles with up to ??l = 2.5 • 10 5 
and N = 4- 10 6 . Details of the algorithmic improvements 
used for this work will be given elsewhere |27j . 



u 

^ 0.5 







— 1 numerical values 

Bulgac et al. 
■ PFA 



le-04 0.001 0.01 0.1 

a/R 



10 



100 



FIG. 2: Casimir interaction energy of a sphere with radius 
R and an infinite plate vs. the curvature parameter a/R. 
The energy is normalized to the zeroth-order PFA formula 
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For larger curvature parameter, the PFA estimate 



(dot-dashed line) differs qualitatively from the worldline re- 
sult (crosses with error bars). Here, we observe good agree- 
ment of our result with the exact solution of |23 [ which is 
available for a/R > 0.1 (dashed line). 



Further advanced field-theoretic methods have been 
developed for Casimir calculations during the past years. 
Significant improvements compared to the PFA have 
been achieved by the semiclassical approximation |28j |. 
a functional-integral approach using boundary auxiliary 
fields [2£j, and the optical approximation [23. These 
methods are especially useful for analyzing particular ge- 
ometries by purely or partly analytical means; in the gen- 
eral case, approximations are often necessary but difficult 
to control. Hence, our results can also shed light on the 
quality of such approximations. 

Sphere above plate. - We consider a sphere with radius 
R above an infinite plate at a (minimal) distance a. A 
contour plot of the energy density along a radial plane 
obtained by a pointwise evaluation of Eq. @ is shown in 
Fig. ^ This density is related to the density of worldlines 
with a given center-of-mass that intersect both surfaces. 
Figure |21 presents a global view on the Casimir interac- 
tion energy for a wide range of the curvature parameter 
a/R; the energy is normalized to the zeroth order of the 
PFA formula ©, E^ A . For small a/R ("large spheres"), 
our worldline result (crosses with error bars) and the full 
sphere- and plate-based PFA estimates (dashed-dotted 
lines) show reasonable agreement, settling at the zeroth- 
order PFA Epp A . The first field-theoretic confirmation 
of this result has been obtained within the semi-classical 
approximation in [2^ ■ The full PFA departs on the per- 
cent level from Epp A for a/R > 0.01, exhibiting a relative 
energy decrease. By contrast, our worldline result first 

stays close to E p °p A and then increases towards larger en- 
ergy values relative to Ej^L. This observation confirms 
earlier worldline studies |22| and agrees with the optical 
approximation pEj in this curvature regime. 

For larger curvature a/R > 0.1 ("smaller spheres"), 
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FIG. 3: Magnified view of Fig. for small a/R. The 
0.1% validity range of the PFA is characterized by curva- 
ture parameters, where the error band of our worldline re- 
sults and the PFA band (blue-shaded/in between the dot- 
dashed lines) overlap, see Eq. JjjJ. The dashed lines depict 
a constraint polynomial fit of the worldline result, p(a/R) — 
1 + 0.35(a/R) - 1.92(a/R) 2 , and its standard deviation. The 
inlay displays the same curves with a linear a/R axis 



we observe a strong increase relative to E p °p A [2(| • Here, 
our data satisfactorily agrees with the exact solution 
found recently for this regime [2j| (dashed line). The 
latter work also provides for an exact asymptotic limit 
for a/R — > oo, resulting in 180/7T 4 for our normalization. 
Our worldline data confirms this limit in Fig. [21 

Two important lessons can be learned from this plot: 
first, the PFA already fails to predict the correct sign 
of the curvature effects beyond zeroth order, see also 
|3C| . Second, the relation between the Casimir effect for 
Dirichlet scalars and that for the EM field is strongly ge- 
ometry dependent. For the parallel-plate case, Casimir 
forces only differ by the number of degrees of freedom, 
cf. the coefficient epp in Eq. Q). For large curvature, the 
Casimir energy for the Dirichlet scalar scales with a~ 2 , 
whereas that for the EM field obeys the Casimir-Poldcr 
law ~ a -4 [32l l3^|. We emphasize that this difference 
does not affect our conclusions about the validity limits 
of the PFA, because the PFA makes no reference to the 
fluctuating field other than the coefficient cpp. 

For a quantitative determination of the PFA validity 
limits, Fig. [3] displays the zeroth-order normalized en- 
ergy for small curvature parameter a/R. Here, our re- 
sult has an accuracy of 0.1% (jack-knife analysis). The 
error is dominated by the Monte Carlo sampling and the 
ordinary-integration accuracy; the error from the world- 
line discretization is found negligible in this regime, im- 
plying a sufficient proximity to the continuum limit. 

In addition to our numerical error band, we consider 
the region between the sphere- and the plate-based PFA 
as the PFA error band. We identify the 0.1% accuracy 
limit of the PFA with the curvature parameter a/i?|o.i% 
where the two bands do no longer overlap. We obtain 



R 



PFA 



0.1% 



< 0.00073 



(5) 



as the corresponding validity range for the curvature pa- 
rameter. For instance, for a typical sphere with R = 
200/^m and an experimental accuracy goal of 0.1%, the 
PFA should not be used for a > 150nm. We conclude 
that the PFA should be dropped from the analysis of 
future experiments. 

For the 1% accuracy limit of the PFA, we increase the 
band of our worldline estimate by this size and again 
determine the curvature parameter for which there is no 
intersection with the PFA band anymore. We obtain 



R 



PFA 



1% 



< 0.00755. 



(6) 



For a sphere with R = 200/im and an experimental ac- 
curacy goal of 1%, the PFA holds for a < 1.5^m. This 
result confirms the use of the PFA for the data analysis 
of the corresponding experiments performed so far. 

In order to study the asymptotic expansion of the nor- 
malized energy, we fit our data to a second-order poly- 
nomial for a/R < 0.1 and include the exactly known 
result for a/R — ► 0. We obtain p(x) ~ 1 + 0.35a; — 
1.92a: 2 ± 0.19x\/l - 137.2a + 5125a: 2 , where x = a/R, 
The fit result is plotted in Fig. (dashed lines), which 
illustrates that E ~ E p °p A p(a/R) is a satisfactory ap- 
proximation to the Casimir energy for a/R < 0.1, re- 
placing the PFA JJJ . The inlay in this figure displays the 
same curves with a linear a/R axis, illustrating that the 
lowest-order curvature effect is linear in a/R. Given the 
results of the PFA QJ, the semiclassical approximation 
pfij . PsAx) — 1 — 0.17a;, cf. j^, and the optical approx- 
imation [2f|, Popt(x) — 1 + 0.05a, the latter appears to 
estimate curvature effects more appropriately. 
Cylinder above plate. - The cylinder-plate configuration 
is a promising tool for high-precision experiments |3ll |. 
since the force signal increases linearly with the cylinder 
length. Figure 2] shows the corresponding Casimir inter- 
action energy versus the curvature parameter. The en- 
ergy axis is again normalized to the zeroth-order PFA re- 



sult, £pp A (a, R) 



-eppjj^ fsTT- The qualitative con- 
clusions for the validity of the PFA are similar to that 
for the sphere above a plate: beyond leading order, the 
PFA even predicts the wrong sign of the curvature ef- 
fects. Quantitatively, the PFA validity limits are a factor 
~ 3 larger than Eqs. owing to the absence of 

curvature along the cylinder axis. 

The most important difference to the sphere-plate case 
arises for large a/R, Here, the data is compatible with 
a log-like increase relative to E p °^ A , implying a surpris- 
ingly weak decrease of the Casimir force for large curva- 
ture a/R — ► 00. Our result agrees nicely with the very 
recent exact result [24| which is available for a/R > 0.1. 
The data thus confirms the observation of [24| that the 
resulting Casimir force has the weakest possible decay, 
F ~ l/[a 3 hi(a/-R)], for asymptotically large curvature 
parameter a/R — > 00. 

In summary, we have computed Casimir interaction 
energies for the sphere-plate and cylinder-plate configu- 
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FIG. 4: Casimir interaction energy (normalized to £"pp A ) of 
an infinitely long cylinder with radius R at a distance a above 
an infinite plate vs. the curvature parameter a/R. The inlay 
shows a magnified view for small values of a/R. 



ration with Dirichlet boundary conditions from first prin- 
ciples for a wide range of curvature parameters a/R. In 



general, we observe that curvature effects and geometry 
dependencies are intriguingly rich, implying that naive 
estimates can easily be misguiding. In particular, predic- 
tions based on the PFA are only reliable in the asymp- 
totic no-curvature limit. Its quantitative validity bounds 
given above and thus genuine Casimir curvature effects 
are in reach of currently planned experiments. 

Beyond the Dirichlet scalar investigated here, it is well 
possible, e.g., for the EM field, that some cancellation of 
curvature effects occurs between modes obeying different 
boundary conditions. In fact, such a partial cancellation 
between TE and TM modes of the separable cylinder- 
plate geometry can be observed in the recent exact result 
for the EM field for small curvature [24| ■ Casimir calcula- 
tions for the EM field in non-separable geometries, such 
as the important sphere-plate case, therefore remain a 
prominent open problem. 

The authors arc grateful to T. Emig, R.L. Jaffc, 
A. Scardicchio, and A. Wirzba for useful discussions. The 
authors acknowledge support by the DFG under contract 
Gi 328/1-3 (Emmy-Noether program) and Gi 328/3-2. 
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